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We present results of a search for the decays f?" (. 7 {H. = e, ji). The search is performed 

using 320 x 10*^ BB pairs collected at the r(45) resonance with the BABAR detector at the PEP-II B 
Factory at SLAC. We find no significant signal and set the following branching fraction upper limits 
at the 90% confidence level: B(B° e+e-7) < 1.2 x 10"^ and ^ IJt^H'l) < 1-5 x 10"^. 

PACS numbers: 13.20.He, 14.40.Nd 



Radiative leptonic decays of neutral B mesons, 
f+f~7 with £ = e, /i [1], are flavor- changing neutral- 
current transitions that are forbidden at the tree level 
in the standard model (SM). In the SM, such processes 
are described by penguin and box diagrams to leading 
order, as shown in Fig. 1. The expected branching frac- 
tions for these processes are of order of 10~-^° [2]. A 
measured branching fraction bigger than this would be 
an evidence for new physics. The largest contributions 
arise from processes in which a photon is emitted from 
one of the initial quarks, thus avoiding the hclicity sup- 
pression of the purely leptonic decay 5° £^£^. A 
search for the processes 5° — > has been performed 

by the BaBAR collaboration and others [3], but there is 
no previous search for the S° i^i^'f decays. 

The analysis described in this Letter uses a sample of 
320 X 10^ BB pairs recorded with the BaBar detector at 
the PEP-II asymmetric energy e"'"e~ storage rings. This 
corresponds to an integrated luminosity of 292 fb~^ col- 
lected at the T(45) resonance. 

A detailed description of the BaBAR detector can be 
found elsewhere [4]. Charged-particle trajectories are 
measured by a five-layer silicon vertex tracker and a 40- 
layer drift chamber operating in a 1.5 T magnetic field. 
A detector of internally reflected Cherenkov light is used 
for charged hadron identification. Surrounding this is 
an electromagnetic calorimeter (EMC) consisting of 6580 
CsI(Tl) crystals, and the instrumented flux return for 
the solenoid, which consists of layers of steel interspersed 
with resistive plate chambers or limited streamer tubes. 

A full BaBAR Monte Carlo (MC) simulation using 
GEANT4 [5] is used to evaluate signal efficiencies and to 
identify and study background sources. The signal MC 
sample is based on a calculation where the B° £+£"7 
transition depends on three Wilson coefficients C7, Cg, 
and Cio at leading order [6]. 

We reconstruct the B° signal candidates by combining 
two oppositely- charged leptons and a photon. The B^ 
vertex is fitted using a Kalman Filter method [7]. The 
leptons are required to originate from a common vertex, 
and the B° candidate is required to be consistent with 
coming from the beam interaction point. 

To minimize the number of misidentified particles, the 
leptons are required to satisfy stringent particle identifi- 
cation criteria [8]. For the electron candidates, the en- 
ergy loss due to bremsstrahlung is recovered whenever 
possible, by looking for the energy deposits (clusters) in 
the EMC close to the intersection of their tracks with the 
EMC. For photon clusters, the transverse shower shape is 



required to be consistent with an electromagnetic shower. 
Leptons and photons are required to reside fully in the 
geometric acceptance of the detector. 

Since the signal event contains two neutral B mesons 
and no additional particles, the total energy of each B 
meson in the ccnter-of-mass (CM) frame must be equal 
to half of the total beam energy in the CM frame. 
We d efine mes = ^/CCJ^^TE^P*F and AE = 
Hi ^^1 + (p*)2 - El^^^, where E^^.^^^ is the beam en- 
ergy in the CM frame, p* and Wj are the momenta in 
the CM frame and the masses of the daughter particles 
i [i = £"'",£",7), respectively, i/boam used instead 
of the measured B meson energy in the CM frame be- 
cause i^beam morc precisely known. For correctly re- 
constructed B" mesons, the rngs distribution has a max- 
imum at the B° mass with a standard deviation of about 
3 MeV/c^ and the /S.E distribution has a maximum near 
zero with a standard deviation of about 30 MeV. 

The B° £+£~7 candidates are selected by re- 
quiring -0.5 < Ai; < 0.5 GeV and 5.0 < rriEs < 
5.3GeV/c2. These ranges include both background- and 
signal- dominated regions. As shown in Fig. 2, five 
background-dominated regions (sideband areas) are used 
for the background estimation. To avoid experimenter's 
bias, the events in the signal-dominated region (signal 
box) and in the shaded area covering the signal box are 
not included in the analysis until the final selection cri- 
teria have been optimized and the background estima- 
tion has been finalized. The shapes of the toes and 
Ai? distributions of the signal MC are parameterized by 
the Crystal Ball function [9] to allow for the asymmet- 
ric shape of the signal peak due to energy loss in the 
EMC. The size of the signal box is chosen to be approxi- 
mately ±3xFWHM for A£; and mEs: -0.146(-0.112) < 
IS.E < 0.082 GeV for the e+e-7 ilJ''^t^~l) mode, and 
5.270 < TOES < 5.289 GeV/c2 for both modes. 

The dominant backgrounds are: 1) unmodeled higher- 
order QED and hadronic two-photon processes for the 
e+e-7 mode; 2) B decays where the photon comes from 
a 7r° decay, or the lepton is from a J/ip or tp{2S) decay; 
and 3) continuum background from e+e- — > // (where 
/ = u, d, s, c, or r) processes at the parton level. 

To take into account higher-order QED and hadronic 
two-photon processes, we introduce additional selection 
criteria for the e+e-7 candidates: we require the cosine of 
the polar angle of e~ (7) to be between —0.743 (—0.618) 
and 0.81 (0.8), the energy of the photon to be > 0.3 GeV, 
the number of charged tracks (EMC clusters) in the event 
to be > 5 (10), and the ratio of the second-to-zeroth order 
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FIG. 1: The penguin (left and middle) and box (right) Feynman diagrams for £'^£~j {£ = e, /i) decays. The photon can 

be emitted from any of the quarks or leptons, but the amplitudes are largest if the photon is emitted from one of the initial 
quarks. 
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FIG. 2: Definitions of the signal box, blinded area (equal 
to the sum of the signal box and of the shaded area), and 
sideband areas in the AE vs. to-es plot: Uppcrl (Ul), Upper2 
(U2), Lowerl (LI), Lower2 (L2), and Middlel (Ml). The 
signal box has the same AE range as the Ml box (different 
for each mode), and the same hies limits as the U2 and L2 
boxes. The figure is not drawn to scale. 



Fox- Wolfram moments {R2) [10], which is calculated with 
the charged tracks and neutral clusters in the rest of the 

event (ROE), to be < 0.7. 

To reduce the number of events where the photon is 
from a tt*^ decay, we veto photon candidates that can be 
combined with any other photon in the event to form a tt*^ 
candidate with a mass within three standard deviations 
(~20MeV/c^) of the nominal 7r° mass. 

We veto Icpton candidates that form a suitable J/^ or 
ipi2S), as described in Rcf. [11]. 

To suppress the continuum background, we require R2, 
calculated from all charged tracks and neutral clusters, 
to be less than 0.35, and the absolute value of the cosine 
of the angle between the thrust axis of the candidate 
and that of the ROE to be less than 0.8. These variables 
arc used in a neural network combined with the following 
variables: 1) the absolute value of the cosine of the angle 
between the direction and the beam axis, 2) the ab- 
solute value of the cosine of the angle between the thrust 
axis of the B^ candidate's decay products and the beam 



axis, 3) the ratio of second order to zeroth order Legen- 
dre moments of all charged tracks and neutral clusters, 
and 4) the invariant mass of the dileptons. The neural 
network rejects 20(36)% of the background while keeping 
95(89)% of the signal, for the e~^e~j {jX^ iJi~^) mode. All 
the selection criteria are optimized with MC samples to 
discriminate signal from background. 

After all requirements are applied, there are on average 
1.01(1.07) candidates per event for the e+e~7 (/i+/U~7) 
mode. In events with multiple candidates, the one with 
the highest probability for the vertex fit is retained. The 
signal efficiency is 7.4(5.2)% for the e+e~7 (/i+/U~7) 
mode. The e"'"e~7 mode has higher efficiency because 
electrons have higher detection efficiency than muons. 

To assess possible background contributions that peak 
in the signal box, we examined 32 exclusive hadronic and 
semileptonic B decays using MC, including events where 
both B mesons decay semileptonically, and found no sig- 
nificant contribution. 

A variety of methods to estimate the background in the 
signal box have been tried, including fitting and counting 
methods in various jties and AE sideband areas with 
different conditions. All studies yield results that are 
compatible within uncertainties. 

The chosen method is model-independent, is based on 
data only, and has a small systematic uncertainty. To 
estimate the background level in the signal box, five dif- 
ferent sideband areas are used, as indicated in Fig. 2. 
The ratio iJ^j is the estimated ratio of the yield in the 
signal box to the yield in the Ml box. The expected 
background in the signal box {'i^g^) is calculated by 
multiplying R^l^ by the yield in the Ml box. We esti- 
mate Rpli as the mean of two ratios R^^ and i?^, where 

box X . This assumes that the changes in the ratio i?^, 
i?^Yt aiid are linear in AE. 

To test our assumption of this linearity, we use MC 
samples and calculate the ratio R^ by dividing the yield 
in the signal box by the yield in the Ml box. The relative 
difference between R^ and R^g^ in MC samples is as- 
signed as a systematic uncertainty. The estimated back- 
ground is 1.75 ± 1.38 ± 0.36 (2.66 ± 1.40 ± 1.58) events 
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for the e+e~7 [ji^ mode, where the stated errors 
represent the statistical and systematic uncertainties, re- 
spectively. 

The dominant source of systematic uncertainty on the 
signal yield is the calculation used for the signal MC [6] . 
The three theoretical input parameters, the Wilson coef- 
ficients C7, Cg, and Cio, used in the calculation are varied 
by ±10%, as recommended by the authors of [6]. This 
variation changes the kinematics of the signal events and 
can thereby impact the detection efficiency. The largest 
relative change in signal efficiency by this variation is 
assigned as a systematic uncertainty. 

We have studied e+e" — > /i"'"/i~7 decays in data to as- 
sess the systematic uncertainty in photon reconstruction. 

The systematic uncertainty from the lepton identifica- 
tion has been determined using an independent control 
sample of J/^f) decays. 

The uncertainty on the number of BB events is 
1.1% [12]. 

The systematic uncertainty related to an imperfect de- 
tector simulation is studied using a control sample of 
B^ ~> J/'ipK'^ events. The same continuum background 
suppression requirements are applied on this sample and 
the signal efficiency is calculated. The relative difference 
in the signal efhciencies between data and MC samples 
is assigned as a systematic uncertainty. 

The systematic uncertainty related to the tracking ef- 
ficiency is determined from e+e~ t^t~ interactions, 
with one r decaying leptonically and the other to three 
charged hadrons. All the contributions to the systematic 
uncertainties are added in quadrature and summarized 
in Table I. 

TABLE I: Summary of the systematic uncertainties on the 
signal yields. 



e 


^e-7 (%) (%) 


Signal Calculation 


2.3 


3.8 


Photon Reconstruction 


1.6 


1.6 


Lepton Identification 


0.7 


1.3 


Number of BB Pairs 


1.1 


1.1 


Data/MC comparison 


1.3 


0.4 


Tracking Efficiency 


0.9 


0.9 


Total 


3.5 


4.6 



After applying the selection criteria we find one event 
in the signal box for each mode, as shown in Fig. 3 and 
Table II. These numbers are compatible with the ex- 
pected backgroimd for both modes. 

An upper limit on the branching fraction is computed 
from 



BuLiB'' ^ £+r7) 



UL 



Nbo ■ ts 



(1) 



the one observed event in the signal box and the es- 
timated background, using the frequentist method de- 
scribed in Ref. [13] including both statistical and sys- 
tematic uncertainties, TV^o is the number of neutral B 
mesons and e^ig is the signal reconstruction efficiency. 
The systematic uncertainties are included in e^ig. It is 
assumed that S(r(4S') ^ B^W) = ^(7(45') ^ B+B"), 
and so is equal to the number of BB events. The 
90% C.L. branching fraction upper limits obtained are 
B{B^ e+e-7) < 1.2 X 10"^ and B{B^ IJ-'^fJ-'l) < 
1.5 X lO-'^. 



TABLE II: Summary of the results where Hobs (nlg^) is the 
number of observed (expected background) events in the sig- 
nal box, €sig is the efficiency, Nul is the 90% C.L. upper limit 

for the signal yield, and Bul is the upper limit on the branch- 
ing fraction at the 90% C.L. The stated uncertainties on n^^^ 
are statistical and systematic, and the uncertainty on €sig IS 
systematic. 



Mode 




nl7 e^rg (%) 


Nul 


Bul 


e"'"e~7 
/^+/^"7 


1 
1 


1.75 ±1.38 ±0.36 7.4 ±0.3 
2.66 ± 1.40 ± 1.58 5.2 ±0.2 


2.82 
2.55 


1.2 X 10"' 
1.5 X lO"'' 



|0.4- 



<0.2 


-0.2 

-0.4 

5" 



5.1 



5.2 
rriEs (GeV/C 




5.2 5.3 
rriEs (GeV/c^; 



where Nul is the 90% confidence level (C.L.) upper limit 
for the signal yield, determined by taking into account 



FIG. 3: Distribution of events in niES and /S.E. Tiie left plot is 
for the e^e^7 mode and the right plot is for the /i^/i^7 mode. 
The dots are the events outside the signal box (rectangular 
region) , and the triangles are the events inside the signal box. 



In summary, a search for B^ — > ^+£"7 (£ = e or 
jj) decays has been performed based on 320 x 10^ BB 
events. We obtain 90% C.L. upper limits for the branch- 
ing fractions of B{B° e+e-7) < 1-2 x lO"'' and 
B{B^ l.i~^fi~j) < 1.5 X lO"*", which represent the first 
limits placed on these decay channels. These are well 
above the SM expectations. 
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